Giant slip lengths of a simple fluid at vibrating solid interfaces 
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It has been shown recently [PRL 102, 254503 (2009)] that in the plane-plane configuration a 
mechanical resonator vibrating close to a rigid wall in a simple fluid can be overdamped to a frozen 
regime. Here, by solving analytically the Navier Stokes equations with partial slip boundary condi- 
tions at the solid fluid interface, we develop a theoretical approach justifying and extending these 
earlier findings. We show in particular that in the perfect slip regime the above mentioned results are, 
in the plane-plane configuration, very general and robust with respect to lever geometry considera- 
tions. We compare the results with those obtained previously for the sphere moving perpendicularly 
and close to a plane in a simple fluid and discuss in more details the differences concerning the de- 
pendence of the friction forces with the gap distance separating the moving object (i.e., plane or 
sphere) from the fixed plane. Finally, we show that the submicron fiuidic effect reported in the 
reference above, and discussed further in the present work, can have dramatic implications in the 
design of nano-electromechanical systems (NEMS). 

PACS numbers: 47.61. Fg, 47.15. Rq, 85.85.-fj, 07.79. Lh 



I. INTRODUCTION 

Nanomcchanical resonators and nanoelectromechani- 
cal systems are widely used in a multitude of applications 
such as ultrafast actuation and sensing at the zeptogram 
and sub-attonewton scales [345| • The extraordinary sen- 
sitivity that such NEMS provide rely mainly upon their 
very high oscillating pulsations and quality factors Q. 
However, while Q factors in the tens or hundreds of thou- 
sands can be obtained in vacuum and/or cryogenic envi- 
ronment, these values degrade dramatically in liquid and 
gas phases meaning that much more work is still to be 
done to reach the technological level (see however 
It is therefore necessary to characterize more precisely 
the viscous forces exerted by fluids on the vibrating mo- 
tion of NEMS and this constitutes the motivation for the 
present work. 

Over the last years, micro- and nano-fluidics experi- 
ments [3] reported that the physical properties of fluids 
flowing into or around confined systems, such as nano- 
channels [sj, are strongly modified compared with those 
encountered at the micro and macro scales. In particular, 
it has been shown that the usual no-slip boundary condi- 
tions, which are universally used since the 19*'' century to 
model the behavior of Newtonian fluids at a solid inter- 
face, break down at the nano scale Such modifications 
of boundary conditions are also expected to have a huge 
impact on NEMS dynamics since properties known for 
Microelectromechanical systems (MEMS) (Tol - [l5| cannot 
simply be scaled down to the nano realm. 
In a recent work, we investigated the importance of gas 
damping on a thermally actuated microlever as it is grad- 
ually approached towards an infinite wall in parallel ge- 
ometry [l6| . The experiment performed at room temper- 
ature in air showed that the sub- Angstrom lever oscilla- 
tion amplitude, i.e., in the direction perpendicular to the 
parallel planes, is completely frozen as the gap d is pro- 
gressively decreased from 20 /im to 400 nm. Moreover, 



the friction force recorded was much larger that predicted 
by the Navier-Stokes hydrodynamical equations solved 
together with the no-slip boundary conditions. Instead, 
the reported results are qualitatively and quantitatively 
well understood if one accepts the perfect-slip boundary 
conditions for which friction at the lever/gas interface is 
prohibited. 

The aim of the present article is twofold. First, we study 
theoretically the motion of the Newtonian fluid, i.e. air, 
around the oscillating micro plate. Starting with the lin- 
earized Navier-Stokes equation we analytically show that 
perfect slip boundary conditions lead to the correct dy- 
namical behavior reported experimentally in ref. [l6j . We 
compare our findings with the other more traditional ap- 
proaches based on the no-slip boundary conditions and 
show that they necessarily conflict with the experimental 
facts. During the analysis we also discuss the different 
possible boundary conditions and in particular the im- 
pact of the slip length of the beam dynamical behavior. 
The second goal of this paper is to show the important 
implications that our findings may have on the engineer- 
ing and architecture of future NEMS operating in gaseous 
environment. Here, on the basis of the results obtained 
in ref. J^] for a model system we discuss precisely the 
existence of a critical overdamped regime for NEMS os- 
cillating in fluids and study the influence of materials, 
geometrical, and intrinsic dynamical properties on this 
regime. 

The paper is organized as follows. In Sec. II we dis- 
cuss the general characteristics of the dynamic problem 
in conjunction with Navier-Stokes equations and bound- 
ary conditions. In Sec. HI the static regime correspond- 
ing to various slip lengths and valid for small gap val- 
ues is studied analytically and compared (successfully) 
with the experimental results reported in ref. [T6j . In 
sec. IV we briefly discuss the possible microscopic mech- 
anisms involved in order to explain the reported results. 
In particular, we compare the calculations obtained here 
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for the plane-plane geometry to those already obtained 
in the s phe re-plane configuration by Taylor and Vino- 
gradova Finally, in Sec. V we discuss the conse- 

quences of our findings for NEMS dynamics. A summary 
is given in Sec. VI. 



II. MECHANICAL OSCILLATOR IN A 
NEWTONIAN FLUID 

We consider a thin silicon commercial cantilever used 
in atomic force microscopy ( AFM) for liquid imaging [l6| . 
This system is modelled as a parallelepiped with length 
L, width w and thickness t. The lever is clamped by 
one end and oscillates mainly along the z direction (see 
Fig. 1). In vacuum this beam may be viewed as a har- 




FIG. 1: Scheme of the one end-clamped parallelepiped can- 
tilever used in ref. [T6[. The lever has dimensions t, L, and 
w and is vibrating in the z direction at a distance d from the 
substrate. 

monic oscillator with an intrinsic resonance pulsation ujq 
and effective mass m (internal losses can be fairly ne- 
glected in the following analysis). In the fluid, the viscous 
force acting on the lever is characterized by a dissipative 
coefhcient 7 connecting the viscous force Fz normal to 
the lever to the velocity U of the lever along the same 
direction: = —7 • U. At short distance, i.e, in the non 
retarded regime, 7 is given by 16]: 



7 



2rjLw 



(1) 



where rj is the (dynamic) fluid viscosity and d the dis- 
tance between the lever and the substrate (see Fig. 1). 
As we will now show this law results directly from the 
perfect slip boundary conditions applied to the Navier- 
Stokes equation, and it differs from the usual 7 ~ 
rjwL^ /d^ behavior deduced by considering the no slip 
conditions [l6|, Il9| . 

To derive this result we start from the non-linear Navier- 
Stokes equations for an incompressible fluid 



dtv + {v ■ V)t/ 



(2) 



where v is the local fluid velocity, p its density, and p 
the gas pressure. To calculate the gas flow around the 



lever we take into account specifical properties of the sys- 
tem under study. First, the importance of the non-linear 
term {v ■ W)v can be estimated from the knowledge of 
the Reynolds number := v ■ Xjv — pv ■ X/rj calcu- 
lated for a characteristic length X and fluid velocity v 
[v — rj/ p — 1.5 X 10^^ m^/s is the kinematic viscosity). 
In the present context, a correct order of magnitude of 
the velocity in the fluid is given by the lever velocity com- 
ponents Ux, Uz '■— U along the x and z direction respec- 
tively. An important related feature is that we have here 
Ux ^ U. Indeed, writing 9 ~ 5z/L (with 6z = 0.05 nm 
the typical lever oscillation amplitude) the main angle 
shown by the lever with the x axis we get the estima- 
tion Ux/U ~ 9 ~ 10^^ which implies that the motion is 
mainly vertical. To simplify our analysis we therefore as- 
sume that the lever is a horizontal plate vibrating along 
the z direction. The lever oscillating at the frequency 
a;/(27r) = 50 kHz one obtains U ~ w6z /{27r) ~ 2.5 x 10"^ 
m/s and thus with X := d = 50 /um 



Re 



■jSz ■ d 



0.8 X lO-'^ < 1. 



(3) 



As a consequence of this ultra-small Reynolds number 
we will completely neglect nonlinearity in the rest of this 
work. As a corollary of this analysis we also deduce the 
Mach number M = u/c ^ J7/c ~ 10^^ ^ 1 (c sound 
velocity in air). The vanishing value of M justifies the 
fluid incompressibility hypothesis V • z7 = 0. 
The second question that we should deal with concerns 
the amplitude of the dynamical term pdtv = —iujQpv in 
Eq. [21 Since we are concerned with harmonic oscillations 
we introduce a second Reynolds number: 



Rl — 



Wo 



2 X IQ-^ ■ (d[/im])^ 



(4) 



where dlfxm] is the measure of d in micrometers. Clearly 
i?g (and therefore pdtv) is negligible as far as d <C 
y^iT/oJo — 10 /um, i.e. as far as the gap d is smaller 
than the boundary layer thickness 5b (see below for dis- 
cussion). In the present work we will only consider this 
static regime and neglect the dynamic term. 
In order to solve the linearized Navier- Stokes equation 



T]W^v~Wp ~ 0, 



(5) 



one must provide the precise boundary conditions for 
the fluid velocity at the solid interfaces. The condition 
for the normal component of the velocity is quite obvi- 
ous d^HO]. Indeed, since the fluid cannot go through a 
solid interface the fluid velocity Vz must equal the veloc- 
ity of the plate U at z = d and must also vanish along 
the surface z = [l^, [lO]. However, the precise form of 
the conditions for the tangential components Vx and Vy 
is not so natural and is a subject of debates and contro- 
versies since the birth of hydrodynamics [2l| - [23l . The 
problem was already well addressed by Navier |24| who 
introduced the two most known possibilities which are 
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respectively the no-slip and perfect slip boundary condi- 
tions. Following the no-slip hypothesis the fluid velocity 
along the interface must equal the in-plane velocity of 
the solid boundary. In the present case the no-slip con- 
ditions lead to Vx ~ Vy — along the substrate plane 
z — and Vx — Ux — 0, and Vy — along the cantilever 
interface at z = d. The no-slip hypothesis is well docu- 
mented in the literature and experimentally justified at 
the macroscale [1, [2l| - |23| . It leads however to increasing 
difficulties and contradictions in the micro and nanoflu- 
idic regime where fundamental and statistical properties 
of molecules such as the mean free path and the surface 
roughness cannot be ignored 0, H, [2l| - [23| . 
The second extreme possibility, the perfect slip bound- 
ary conditions, suppose having dzVi = (from now on we 
use the notation i = x,y) along the solid interface. These 
are reminiscent for conditions on the viscous stress tensor 
= r]{dzVi + diVz) — meaning that no tangential fric- 
tion is allowed between fluid and solid. We however point 
out that the equivalence between the relations <7iz = 
and dzVi = assume the additional conditions diVz — 
which implies that Vz does not depend on x, y in the 
vicinity of the solid boundary (this is obviously true on 
the interface provided that the boundary conditions on 
Vz are fulfilled). Beyond these two extreme cases more 
realistic approaches were proposed to take into account 
a possiblepartial slip at the boundary. In particular 
Navier 0, H, l2l| - |23 | suggested the existence of a surface 
friction force law Uiznz = nvi where fi = UzZ is a, unit 
vector normal to the surface and oriented outwardly from 
the solid to the fluid and k a dissipative coefficient. With 
the same assumptions as before this leads to 



rizdzVi = KV^/ri = vjb, 



(6) 



where b is the so called slip length. This law has been 
considerably studied in the recent years with the ad- 
vent of micro and nanofluidics 0, M, HI, HSl- Eq. [B] 
has been confirmed many times in particular for liquid 
flows in nanochannels |22| . However, the measurements 
of the associated slip length b reveals a broad spectrum 
of values which specifically depend on the system consid- 
ered 0, . This shows that only a microscopic analysis 
could lead to a better understanding of the phenomenon. 
Keeping this point for latter discussions we will here ap- 
ply the law given by Eq. |6] to our problem and see how 
it compares with the experimental results. 



III. THE STATIC LIMIT 

In the present analysis we consider the static regime 
valid for d ^ 6b ■= \j2v jiOQ. We must therefore solve 
the system of coupled equations TyV^t? = Vp, V • w = 
together with the conditions given by Eq. [S] The prob- 
lem is reminiscent for the one solved by Reynolds in the 
case of two parallel disks in a dissipative fluid. Reynolds 
considered two disks approaching each other with a con- 
stant velocity ±C/ along their common axis of symmetry. 



However, despite geometry differences, Reynolds consid- 
ered speciflcally the case of no-slip boundary conditions 
(which were universally accepted at that time) and not 
the more general Eq. [B] Using the same approxima- 
tions than Reynolds we here suppose dxVi,dyVi <^ dzVi, 
dxVz,dyVz dzVz, and dzP — which are standard in 
lubrication theory. We then have 

f]d'^Vi ~ dip. (7) 

After integration with respect to z this leads to 

Vi{x,y,z) = ■^dip{x,y)z^ + ai(x,y)z + /3i{x,y). (8) 

To be general we are introducing two a priori different 
slip lengths bo and bi for the interfaces at z = and 
z = d, respectively. This hypothesis implies 



Vi{z 0) = bodzVi{z = 0), Vi{z ^ d) 
and therefore 



-bidzVi{z = d), 



v,{x, y, z) = ^dip{x, y)[z'^ - d ^^^'^^ , (z + &o)]- (9) 

The fluid incompressibility relation tiaj+SyWy-f^^w^ — 
and the boundary conditions for Vz at z — Q give the 
expression: Vz{x, y,z) = - J^{dxVx + dyVy)dz, i.e.. 



vz^^{d'x + dl)p{x,y)[d 



2bi + d ,7? ^ ^ z3 



(10) 



At z = c? we have the boundary condition Vz = U and we 
deduce 



idl + d^)p{x,y) 



2r]U 



and 



Vzix,y,z) 



(11) 



(12) 



In particular, in the limit bo = bi ^ +oo we have 
Vz = Uz/d whereas for bo = bi = we obtain Vz = 
6U{-z^/3 + z^d/2)/d^. It should be observed that the 
solution for Vz in the perfect slip limit looks like the well 
known solution of the Couette problem [l^ IIO] for the 
permanent fluid motion between two plates in relative 
and uniform motion along the x direction. It is worth 
noting that the direction of the fluid motion and the 
boundary conditions used are however completely differ- 
ent in these two problems (indeed, in the Couette prob- 
lem we assume the no-slip boundary conditions and ne- 
glect Vz [iMIi^l)- Here due to the presence of terms with 
Vz our approach deviates from this standard result. 
The normal force exerted by the fluid located between 
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z = and z = d on each surface element dxdy of the 
lever is given by 

(13) 

Gzz = '^'ndzVz is the fluid stress tensor along the z direc- 



tion and therefore — cr. 



idxdy is the dissipative part 



of the resulting force dFz- Using Eq. [T2]we obtain 



=d = 



2bi+d 



1 2fci+d N .3 

2 60+bi+d''" 



2bi+d 
bo+bl+d 



(14) 



Interestingly this contribution is independent oi x,y and 
of the lateral boundary shape (C) associated with the 
lever and the substrate. The resulting dissipative force 
— ^f^g^crzz\z=ddxdy on the plate of surface S is conse- 
quently 



2'qSU[d^ 



2bi+d 

bo+bi+d 



(d^ + bod)] 



1 2bi+d 

2 bo+bi+d 



)d3 



2bi+d 
bn+bi+d 



bod^] 



(15) 



The second contribution to the force is associated with 
the volume pressure p{x, y) that is solution of the 2D 
Poisson equation [d^ + dy)p{x, y) = A where A is a con- 
stant (see Eq. 11). Lateral boundary conditions along 
(C) are here playing explicitly a role in the analysis. Let 
us suppose for example that the substrate is infinitely 
extended in the z = plane and that the lever is delim- 
ited by a rectangular boundary of length L and width 
w in the z = d plane. A lengthy Fourier analysis (see 
appendix) shows that the pressure field can be expanded 



p{x,y) ^po 



AAL^ 



y 4t sin (7rnj;/L)[a„e+'^"^/^ 
z — / 

n 

+/3„e-^"''/^-l],(16) 



where n are odd integers and , /3„ are constant coeffi- 
cients (see appendix). To obtain this formula the bound- 
ary condition p — pq [pq is the atmospheric pressure) was 
imposed on the rectangular contour. Similarly, if we con- 
sider the ID problem with L finite, w = +00, and with p 
a function of x only we deduce after a direct integration 



Pix) = Pa 



A, 2 
2^" 



xL), 



(17) 



where we used the boundary conditions p(0) = 
p{L) = po- Equivalently, this result could also be 
obtained from the Fourier expansion p{x) = Pq — 
X^n s\n{'Knx/ L). Finally, as a last example, we 
consider the circular plate of radius R. Assuming the 
radial symmetry this problem leads to the solution 



p{r)=pQ + ^{r^ -R"), 



(18) 



where r is the radial coordinate and p{R) = Po along the 
circular contour of radius R. More generally, by using the 



uniqueness theorem it is easily shown that the pressure 
p{x, y) can always be written 



P{x,y) ^ pa + A- g{x,y), 



(19) 



where g{x,y) is the solution of {d^ + dy)g{x,y) = 1 
which is univocally determined by the boundary condi- 
tion g{x, y) ^ along the contour. The resulting vertical 
force due to the pressure field on the plate of surface S 
is thus given by 



^pressure 

= A 



(S) 



{p{x,y) -pQ)dxdy 

(S) 

g{x,y)dxdy = AS{g), 



(20) 



where the term —po equilibrating the pressure -|-po from 
Eq. 19 is due to the force exerted by the fluid on the 
second side of the cantilever. Considering the examples 
quoted previously we find explicitly 



rcssurc 



8AL^ 



n 

irnw fi+TTnw / L ^ — TTnw/L ^ 

for the rectangular plate, and 



pprc 



7]-4 „4 



AL^ 



12 



(21) 



(22) 



for the ID plate (note that w is a finite width in the 
y direction introduced for reasons of dimensionality and 
that Eq. 22 is the hmit of Eq. 21 for w — > +00). Similarly, 
for the disk we deduce 



^prcssu 



7rAR4 
8~~ 



(23) 



These forces due to pressure can be compared with the 
dissipative contribution given by Eq. 15 which can equiv- 
alently be expressed as 

j^dissip. ^ ^^[^2 _ ^^1+1^(^2 ^ (24) 

60 -I- Oi -I- d 

and leads to the total force F, ^ ppressmc ^ ^dissip.. 



Fz^AS[{g)+d' 2^^ + ^ -2 



with 



A = 



bo + bi + d 
2riU 



{d' + bod)l (25) 



(26) 



3 ^ 2bo+bi+d)^ ^ fco+bi+rf ° ^ 



In the usual limit of vanishing slip lengths F^issip. _ q 
and we therefore have 



d3 ^ 



(27) 
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which is the generahzation of the Reynolds formula for 
an arbitrary plate shape in the no-slip limit. Inversely, 
in the limit 60 = ^1 infinite, i.e, perfect slip, we obtain 
y)=pn, FP'"^^''''' = and thus 

which is the result used in ref. in the perfect slip limit. 
It is worth noting that this formula is only dependent on 
the surface S of the vibrating plate and not on its exact 
geometry. It implies that the result obtained in ref. [l^ 
should be very robust with respect to geometry consid- 
erations. We point out that this perfect-slip limit can 




— b=+cc 

— b=100t[m 
b=500nm 

• experiment 



d[^m] 



FIG. 2: Evolution of the beam quality factor Q as a function 
of the gap d for different slip lengths 6 and for the lever me- 
chanical properties given in ref. The experimental data 
are in good agreement only with b > 500 fim (red curve). 



From this follows the damping constant 7 := —Fz/U. 
To compare this result with the experiment we show on 
Fig. 2 the evolution of the quality factor Q — fc/(ajo7) of 
the oscillating beam calculated for a stiffness k = 0.0396 
N/m and a pulsation ujq — 2tt x 50 kHz considered in 
ref. [l3] and for different slip lengths b. Clearly, the agree- 
ment with the experimental data is very good in the limit 
d < Sb (i.e, far from the observed saturation at large 
gaps) for giant slip lengths such as 6 > 500 /xm (red 
curve). Oppositely the reported results are in total con- 
flict with the standard no-slip prediction Q ^ (P (blue 
curve). While the present study focussed on the static 
regime, it is interesting to remark from dimensional anal- 
ysis that one should expect to observe saturation around 
5b — 10 /xm (i.e., 7ii,„ = 2rjS/SB)- The value observed 
experimentally corresponds to Sb — 25 /im, which is 
of the same order of magnitude but nevertheless signifi- 
cantly larger. The difference could be imputed to geom- 
etry considerations, i.e., to the fact that the dynamics of 
the fluid around the lever should be strongly influenced 
by the flnite size of the system under study (indeed, in 
the regime where the boundary layer Sb plays explicitly 
a role, retardation should be taken into account). A dif- 
ferent explanation could be that the model of perfect slip 
breaks down at large gap. A preliminary analysis in that 
direction shows that if we conserve the inertial term in 
the Navier-Stokes equations we indeed obtain a satura- 
tion regime due to an additional damping occurring on 
the length scale b ^ Sb and this even if b is inflnite in the 
static regime. 



be directly obtained by solving the equation rjO^Vz = 
with the boundary conditions for v^- Indeed, the direct 
integration gives Vz = Uz /d. The pressure p — po can be 
deduced similarly [25'|. 

It is worth noting that in the derivation of all these for- 
mulas the Reynolds assumption dzP — implies also 
TjdzVz = dzP = 0. However, from Eq. 12 we deduce 



2bi+d d^ 



d2' 



(29) 



2' bo + bi 

In the perfect slip limit (where ^ 0) the Reynolds 
condition is automatically fulfllled, i.e., the solution is 
self-consistent. However, in general rjdzVz ^ and this 
limits the validity of the previous results. The Reynolds 
hypothesis can nevertheless be justified, for all practical 
purposes, if \d^zVz\ « that is, if \z - ^l^^f | « 

\dig\. Since \dig\ ^ L, where L is a typical lateral dimen- 
sion of the lever, the Reynolds assumption is in general 
justified for L ^ d. 

To study the infiuence of the slip length on the dynamic 
properties of the oscillator considered in ref. we will 
now consider a simplification. Such a simplified force 
expression is obtained if the condition 60 = ^1 = ^ is ful- 
filled. Considering for example the rectangular plate (in 
the limit given by Eq. 22) we get for the total force 



IV. DISCUSSION 

Fundamentally, the existence of a giant slip length 
regime is very surprising and interesting and should 
therefore be discussed carefully. Here, we will only review 
some results which, we think, are important to justify 
microscopically the results discussed in this article and 
in ref. [16|. We remind that from a microscopic point of 
view, the slippage coefficient is actually linked to the very 
nature of the interaction between the oscillating surface 
and the air molecules. Historically, the first theoretical 
analysis of this phenomenon goes back to J. C. Maxwell 
and to its kinetic theory of gases [2^. Following this 
approach one can indeed distinguish between a specular 
channel of interaction, for which the molecules are col- 
liding elastically with the surface, and a channel of inter- 
action for which molecules are reflected diffusively from 
the wall [i,[2l[2ii|. This second channel is linked to 
multiple collisions between molecules and also to adsorp- 
tion by the surface. The slip length b in this statistical 
model is given by the Maxwell formula [2l|, [26l - [30| : 



3 Pd 



Pd 



(31) 



where A is the typical mean free path of gas molecules 
(i.e., A ~ 60 nm for air in ambient conditions) and pd the 
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tangential momentum accommodation coefficient, i.e., 
the fraction of tliose molecules hitting the surface which 
are reflected diffusively. Clearly, if pd vanishes then the 
slip length is infinite. This suggests that in the work- 
ing regime of our mechanical oscillator the molecules are 
mainly reflected specularly. Furthermore, recent anal- 
ysis based on the fluctuation dissipation theorem and 
the Green-Kubo relationship emphasize the importance 
of several other microscopical parameters on the molec- 
ular dynamics close to a surface 0, [H, . Altogether 
these studies provide an estimate for the slip length given 
by: 





plane-plane 


sphere-plane 


b = 




7 = QirriR'' /d 


b — ^ +00 


7 = 2riS/d 


7 


= 27rriR^ In {d/{6b))/b 



ksT-qD 



(32) 



where D is the fluid diffusion coefficient, cr, e are re- 
spectively a typical length and energy characterizing 
the molecular interaction, and C_l a coefficient measur- 
ing the roughness (larger C_l mean larger roughness). 
This shows the importance of surface roughness but also 
of surface defects and nano-structuration [ll] for the 
physics of slippage at solid-gas interfaces. Future work 
should investigate the effects of theses parameters on the 
damping coefficient 7. It is worth noting that past studies 
on the slippage at a solid-fluid interface mainly focussed 
on liquids for which the mean free path is much smaller 
(A ~ 1 nm) and for which interactions between molecules 
are much stronger. The existence of a partial slip regime 
implies in those cases to work in the realm of nanoflu- 
idics (e.g., nanochannels) with separating gap d well be- 
low the micrometer range |7bi22,] or with ATMs in contact 
mode [13, HH . Here oppositely, we consider gases and we 
can define a Knudsen coefficient Kn = X/d ~ 0.001 — 0.06 
which corresponds to a regime of transition flow 0, [1, 
occurring at large gap values d (it is indeed well known 
that important deviations to the no shp boundary con- 
ditions appear for iir„ 10^^ - 10"^ Jf). 
Another relevant length in the analysis is the vibration 
amplitude of the lever which we reported in ref. to 
be dz ~ 0.05 nm. This is actually a very small ampli- 
tude which is attainable experimentally mainly because 
of the thermal excitation mode and high sensitivity of 
the optical detection setup used in [l^. This value for 
6z is also comparatively smaller that those attainable 
with actuated AFM [S^l- Working in this regime where 
5z/X ^ <C 1 could therefore lead to new physics 

and we expect that further studies in this direction will 
be done in a close future (e.g., to compare the effect of 
the vibration amplitude and of the excitation modes on 
the micro- and nano-lever dynamics). 
In this context it is worth mentioning that the theo- 
retical model developed in this article for the plane- 
plane configuration constitutes the equivalent of the 
Vinogradova formula obtained for the sphere-plane con- 
figuration i, [H [il m m. The analysis of Vino- 
gradova generalizes the result obtained by Taylor 22 1 
in the no slip limit and which predicts a friction force 
Fz = —GTrrjR^U/d for a sphere of radius R moving along 



TABLE I: Table summarizing the asymptotic viscous force 
regimes for both the plane-plane and sphere-plane configura- 
tion and for the no-slip and perfect slip cases. 



the z axis with the velocity U perpendicularly to the 
interface z = and separated from this plane by a min- 
imal gap d. It is interesting to point out that the no- 
slip condition predicts the same law i^^ cx 1/d as in the 
plane-plane configuration with the perfect slip condition 
(i.e., Eq. 28) but with a different numerical value for 
7Tayior := &TT'qR^ / d (compare Eq. 1). The Vinogradova 
model predicts oppositely ~ 2Trri R^U In {6b /d)/b in 
the limit b +00 of the perfect slip. The different 
regimes of force depending of the value for the slippage 
length b and from the geometry considered are summa- 
rized in Tab. 1. 

The present problem reminds us a very known similar 
difficulty encountered in experiments for measuring the 
Casimir force in the sphere-plane or plane-plane con- 
figuration. It is worth pointing out however, that in 
the Casimir effect [s^l the force in the plane-plane con- 
figuration varies as Fz ~ hS/d'^ whereas it varies as 
Fz ^ hR/d^ for the sphere-plane configuration {h is the 
planck constant). Beside the important difference in the 
power law behavior in d between these expressions and 
those predicted by hydrodynamics it is interesting to ob- 
serve that in the perfect slip limit the viscous force decays 
slowly when d increases whereas the same force vanishes 
in the sphere-plane as predicted by the Vinogradova for- 
mula for b = +00. Therefore, in order to observe the 
perfect slip regime the plane-sphere configuration would 
be much more demanding than the plane-plane config- 
uration studied in this work and in since it implies 
that one should consider very small gaps d to obtain a 
finite effect with the sphere. This could have implica- 
tion in optical near-field microscopy where viscosity is 
suggested as a possible mechanism to justify the shear 
force applied on the tip probes [s^ (see also ]^ and [s^ 
for experimental demonstrations in high vacuum at cryo- 
genic temperatures) . Additionally, this sensitivity of the 
force behavior with b in the sphere-plane configuration 
could be used to probe more precisely the value of the 
slippage length than in the plane-plane geometry. Oppo- 
sitely, the giant slip effect studied in this paper could thus 
be a specificity of the planar geometry and is expected 
therefore to have a huge impact on the NEMS dynamics 
which are mostly developed with such geometry. 



V. IMPLICATIONS FOR NEMS 
ARCHITECTURES AND DYNAMICS 

This brings us to the second point that we shall now 
discuss in this section, that is the implication of our re- 
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suits for NEMS and MEMS engineering. Owing to the 
formula given by Eq. 1 the effect of dissipation on the 
beam motion decays very slowly with d and in particu- 
lar we see that the quality factor Q :— towq/t decreases 
linearly with d. Additionally, when d decreases the pul- 
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10' 



Q. 

8 10' 



" 10-^ 



10" 




10^ 10' 10' 10' 10' 

Resonance frequency [Hz] 



FIG. 3: Evolution of critical gap as a function of frequency 
fo = ajo/(27r) for a Si (black line), SiC (red line), GaAs 
(green line), and Au (blue line) cantilever, respectively. The 
black dashed line corresponds to the limit associated with the 
boundary layer. Experimental data points for the lever stud- 
ied in ref. at the internal resonance frequency 50 kHz are 
also shown for distance d close to the overdamping regime 
(black crosses). Red crosses and circles correspond to char- 
acteristics of NEMS realized with a distance to the substrate 
constantly equal to d = 400nm. Such NEMS can approxi- 
mately be modeled as the one shown in Fig. 1. 

sation at resonance 



(33) 



is progressively down shifted jlQI] . This occurs until 
'^loson. = 0, i.e. from Eqs. 1,33 when the critical distance 
dc given by 



dn = 



V2r]L 



w 



mojo 



(34) 



is reached. For d < dc the lever motion is consequently 
overdamped (nonlinear effect are also expected in this 
limit ). To describe quantitatively the importance of 
such a regime on beam dynamics we remind that for the 

levers considered here we have: ujq — ^ -^fo ' and 
TO — ptLw where p is the bulk density, E is the Young's 
modulus, and vq ~ 1.875104 

Fig. 3 shows the variations of dc as a function of fo — 
ujQ I (27r) for a thin lever with t = 180 nm (same as in 
ref. [l6|) and for different commonly used materials. For 
low frequencies in the 10 kHz range and below, the over- 
damping regime appears already at large separation dis- 
tance dc > 0.1-1 /im. Oppositely, for very high oscillator 



frequencies in the 100 kHz range and beyond we have 
dc ^ 10 — 100 nm and the overdamping regime becomes 
a fundamental issue only at the nanoscale. For compar- 
ison we show on the same graph the experimental data 
points taken from ref. [l6| and corresponding to working 
distances d which are decreasing until the overdamped 
regime at dc is reached. Additionally, we show also the 
physical characteristics recorded (i.e., internal resonance 
frequency /o, and fixed distance gap d = 400 nm with 
the substrate) of typical Si made NEMS. These NEMS 
can be with a good approximation described with the 
simple geometry considered here. Clearly, working with 
such NEMS in a gaseous environment may strongly af- 
fect their dynamics and only for very high /o could the 
overdamping regime actually be overcome. This is in- 
deed confirmed for those NEMS annexed by a red cross 
in Fig. 3 which we studied experimentally by using the 
same optical method as described in ref. . The experi- 
ment showed that NEMS with such gaps do not resonate 
in air at room temperature confirming therefore the role 
played by overdamping. For the NEMS indicated by a 
red circle, i.e., with frequency in the 100 kHz and MHz 
range we were out of the detection sensitivity of our setup 
and no data were available. 
As a complementary analysis we show in Fig. 4 the 
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FIG. 4: Evolution of the critical gap as a function of the ratio 
t/L for the same lever materials as in Fig. 2. Colors lines and 
data points have the same meaning as in Fig. 2. 

explicit dependence of dc on the aspect ratio t/L for dif- 
ferent materials. In the typical range of aspect ratio con- 
sidered the overdamping regime covers distance gaps d 
going from the micrometer range to the nanoscale and, 
therefore, cannot be neglected. Again, this fact is con- 
firmed by comparing these graphs with available experi- 
mental data (see Fig. 4 and compare with Fig. 3). 
Furthermore, it is also useful to remind once again that 
for a large gap d the important length scale is the bound- 
ary layer thickness Sb which characterizes the spatial re- 
gion surrounding the lever for which viscosity has an im- 
pact on the fluid dynamics [l^ [2^. For d > 6b the 
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substratum lies outside this layer and dissipation must 
saturate [l^ [13] , as reported in ref. [11] . Comparisons 
with values for dc (see Fig. 3) show that the overdamping 
regime is always reached for gaps smaller than Sb- The 
overdamping regime appears consequently as a robust 
limitation which should affect the design of any NEMS 
working in fluids. The results obtained here for a par- 
ticular lever geometry are expected to be very general as 
soon as the beam geometrical dimensions are larger than 
6b- However, when dimensions are smaller, boundary ef- 
fects due to the finite size of the system should explicitly 
be taken into account in the definition of dr- 



VI. CONCLUSION 

In this article we studied the linearized Navier-Stokes 
equation in the static regime to describe the damping 
mechanism of oscillating micro plate in air close to a sub- 
strate. We considered the influence of the slip length and 
showed that results reported in ref. Tg'] are only compat- 
ible with very large slip length in the range b > 500/im. 
We discussed the implication of this mechanism on the 
oscillation properties of NEMS and showed that an over- 
damping behavior represents a fundamental mechanism 
for such systems. We expect that this work could have 
important consequences for NEMS engineering. 
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Appendix A: The rectangular plate in the static 
regime 

The aim of this appendix section is to solve the Poisson 
equation {d^ + dy )p{x, y) = A where A is a constant for 



the rectangular domain L x w with boundary conditions 
p = Po along the contour. Using a Fourier series satis- 
fying these boundary conditions at a; = and x — L we 
can write 



p{x, y) ^po + Yl ^"(y) ^™ (" 



(Al) 



Here c„(j/) is solution of 



K-) c„{y) + -T^cn(y) = — e, 
L dy^ rnr 



(A2) 



with e„ = 1 if n is a odd integer and e„ = if n is even. 
The general solution of this equation is 



4A 

-j3/7T-3 



(±) ■ 
The constant an are determined by the boundary con- 
ditions p = Po a,t y = and y = w. One therefore obtain 



(-) 



Jg+irnto/L _ 2J 



This leads to 



r2 



[1 — g-'^"'"j/ij 



(A4) 



p{x,y) =_po + 



AAL 



3 ^^sin(7rna;/i)[-l 



+ 



■K" '■ — ' n 

n 

2^ ^ — -Knw/L 

^+-iTnw/L _ 2 
6~ 



(A5) 
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